A Pulse Shaping Algorithm of a Coherent Matter Wave 
Controlling Reaction Dynamics. 
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A pulse shaping algoritiim for a matter wave witli tlie purpose of controlling a binary reaction has 
been designed. The scheme is illustrated for an Eley-Rideal reaction where an impinging matter-wave 
atom recombines with an adsorbed atom on a metal surface. The wave function of the impinging 
atom is shaped such that the desorbing molecule leaves the surface in a specific vibrational state. 

PACS numbers: 82.53.Kp, 03.75.-b, 03.75.Pp, 32.80Qk, 34.50.Dy 



A binary chemical reaction can be controlled by means 
of interference of matter waves. A source of coherent 
matter waves for such a control scheme could originate 
from an output coupler of a Bose-Einstein condensate^, 
3 0j ■ A two-pulse scenario for such a control has been 
xeviously explored and termed (2PACC) spectroscopy (il 
In the 2PACC model the control knobs were limited 
to the time delay between the pulses and the phase dif- 
ference. R was demonstrated that the control was able 
to obtain a significant enhancement of the probability of 
creating a gas phase molecule in an Eley-Rideal reaction. 
The drawback of the 2PACC scheme was that the desorb- 
ing molecule possessed a broad distribution of vibrational 
states. The present study demonstrates that additional 
control on the incoming matter wave can create a specific 
final vibrational state. 

To achieve this target the task is to figure out how to 
shape the incoming wavepacket. An iterative computa- 
tion scheme is developed for this task which is in analogy 
to the optimal control theory of optical pulse shaping 
[T], I3, y, |lO| . As important is the development of experi- 
mental schemes able to shape such a wavepacket. 

As in earlier studies j^Q, the control target is the Eley- 
Rideal bimolecular reaction where an impinging atom re- 
combines with an adsorbed atom on a metal surface to 
form a diatomic molecule 



Y + H/Cu(lll) ^ YH + Cu(lll) 



(1) 



The source of atom Y is a matter wave, which in the pre- 
sented model consists of either hydrogen or alkali atoms. 
The matter-wave of Y is directed to a Cu(lll)-surface 
with low coverage chemisorbed hydrogen atoms. When 
the wave function of Y overlaps with the one of the ad- 
sorbed atom, interaction is expected, leading to a re- 
combination that forms the YH molecule. H the newly 
formed molecule has sufficient kinetic energy it will even- 
tually desorb from the surface. By shaping the wave func- 
tion, the vibrational state of the desorbing molecules be- 
comes controlled by constructive or destructive quantum 
interferences. 

For a two electronic state case, the wave function of 



the system is described by a vector 



*(r, Z, t) 



z,t) 

*p(r, Z,t) 



(2) 



where the reactant and product channels are denoted by 
the index {i?, P}. The two degrees of freedom are the in- 
tramolecular distance between the two atoms (r) and the 
distance from the surface to the center-of-mass motion of 
the molecule {Z). 

The dynamics of the Eley-Rideal reaction was followed 
by solving the time-dependent two-channel Schrodinger 
equation. The pulse shaping algorithm is based on prop- 
agating the wave function forward (t > 0) and backward 
(r < 0) in time. The Chebychev propagation method[ll| 
was employed since it allows the propagation of the wave 
function in both time directions 



*(r, Z,t + T) 
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afcK)*fc(r,Z). (3) 
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FIG. 1: (Color) An initial trial wave function is shown for 
the pulse shaping scheme for the Eley-Rideal reaction where 
an impinging lithium recombines with chemisorbed hydrogen 
atom on a Cu(III)-surface. This wave function is used for 
shaping the matter wave which leads the molecule to desorbs 
in the fifth vibrational state (i;=4). 



2 




4 6 8 
r in Angstrom 




3 4 5 6 7 
r in Angstrom 



FIG. 2: (Color) Above, the selection step of the pulse shaping scheme is illustrated. A wave function on the reactant channel, 
(/>^, is shown after termination of the backward time propagation (left). The wave packet of the red box has been selected to 
represent the wave packet of the matter wave on the reactant channel. The renormalized selected wave function, which 
represents the wave function of the matter wave is shown (right). This wave packet is used for the backward propagation. 
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FIG. 3: (Color) Above, the projection step of the pulse shaping scheme is illustrated. The trial wave function before the 
projection is showed on the left. The wave function on the product channel, </)5{, after termination of forward time propagation 
(center). The trial wave function after projection (right). 



The coefficients of the expansion, Uk , are proportional to 
the k'th Bessel functions Jk- These coefficients depend 
on the sign of the time propagation. The time symmetry 
enters in the Bessel functions Jk, the Bessel functions of 
odd order change sign as the direction of time propaga- 
tion changes. The functions, are obtained by op- 
erating Tfc(i?jv) on 'S'(r, Z, i), where are the Cheby- 
chev polynomials. The normalized Hamiltonian Hn is a 
shift-and-scaled version of the original Hamiltonian, en- 
forcing the eigenvalues of the Hamiltonian into the inter- 
val [—1, 1]. The parameters, co and i^, are a result of the 
normalization of the Hamiltonian. 

The task is to find the initial wavepacket on the reac- 
tive channel which will lead to desorption of the molecule 
in the w'th vibrational state. The scheme is based on re- 
versing the reaction e.g. starting from a product wave 
function representing the desorbing molecule in the t;'th 
vibrational state. This wave function is propagated back- 



wards in time to find the reactive wave function lead- 
ing to the desired final condition. The scheme is illus- 
trated for an Eley-Rideal reaction in which an imping- 
ing lithium atom recombines with an adsorbed hydrogen 
atom. A trial wave function on the product channel has 
been constructed. Far enough from the surface the in- 
ternal and external degrees of freedom are uncoupled. 
Thus the trial wave packet can be expressed as a prod- 
uct of the w'th molecular vibrational eigenstate, Xvif) 
{v—{0, 1, ...}), and a Gaussian wavefunction, g{Z), 

The wavefunction used for control is the wavefunction 
of a lithium matter wave source. The target is to des- 
orb a lithium hydride molecule in the fifth vibrational 
state (v — A). A target trial wave function has been con- 
structed according to EqO and is shown in Fig. ^ An 
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FIG. 4: (Color) The wave function on the reactant channel which leads to 67 % of molecules desorb in the fifth vibrational 
state (left) Fourier transform of the a one-dimensional wave function represent the impinging atom for fixed distance between 
the surface and the chemisorbed atom at 0.916A (right). 



iterative scheme is used for inversion i.e. finding the con- 
dition for the wavepacket on the reactant channel lead- 
ing to desorption of the molecules in the w'th vibrational 
state. 

1. Backward propagation: A trial wave function 
^triaii'i^, Z) is propagated backward in time. As 
the molecule approaches the surface population is 
transferred to the reactant channel through the 
non-adiabatic coupling. The part of the wave 
packet which is transferred to the reactant channel 
corresponds to a molecular dissociation. That is 
one atom is chemisorbed on surface and the other 
desorbing to the gas phase. After a preassigned 
time the propagation is stopped, leading to the 
wave function = {cf)^j^, cj)p) . Fig. |21(left) shows 
a wave packet in the reactant channel after termi- 
nation of the backward time propagation. 

2. Selection: Part of the backward propagating wave 
function, *i is selected consisting of a wavefunc- 
tion outside the interaction region. In Fig. [3 (left) 
the wave function is shown before selection. The 
part of the wave function which is inside the red box 
has been selected to resemble the wave function of 
the reactive channel. The selected wave function is 
renormalized and it is shown in Fig. [21(right). The 
size and location of the box is chosen such that the 
energy of the selected and normalized wavepacket 
is close to the energy of the trial wavepacket. 

3. Forward propagation: The selected wavefunction, 
*2 = (</'_R,0), is propagated forward in time. 
Through the non-adiabatic coupling between two 
potential energy surfaces. The newly formed 
molecules desorbs leading to the wavefunction 

4. Projection: The overlap between the forward prop- 
agated wave function, '4'3, and the trial wavefunc- 
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FIG. 5: (Color) The accumulated flux (red line) is shown 
as a function of the iteration number. The percentage of the 
accumulated flux in the forth and fifth state have been com- 
puted as a function of the iteration number. The forth and 
fifth vibrational states are shown with triangles and diamond, 
respectively. 

tion is evaluated as a function of the center of mass 

c„(Z) = J ^tr,ai{r,Zy^3ir,Z)dr . (5) 
A new trial wavefunction is created 

^) ^ ( c„(Z)/(Z)xn (r) ) • 

In the first iteration f{Z) is g{Z). This new wave 
function is renormalized and it is used as a new 
trial wave function in step 1 "Atrial — "^proj- The 
projection of the forward propagated wave packet 
is illustrated in Fig. O 

The above pulse shaping scheme is repeated a number 
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of times. It is important that the energy of the total sys- 
tem is conserved after either the selection stage or projec- 
tion of the total wave function. The pulse shaping algo- 
rithm is stopped when a satisfactory yield in the desired 
vibrational state has been reached. The ultimate goal 
is that all the flux of desorbing molecules is probed into 
the particular vibrational state. In order to improve the 
convergence of the algorithm the number of grid points 
has been enlarged (1024x1024) and the grid spacing has 
been reduced {dr—dZ—Q.05a.u.) compared to the earlier 
2PACC calculations |Q] . The potential energy surfaces are 
given in Ref.[^, where the non-adiabatic coupling term 
was increased to /3=0.1eV. The backward propagation is 
stopped when the non-adiabatic population transfer is 
negligible. The size and location of the box used for se- 
lecting the wave packet in step 2 vary from iteration to 
iteration in order to keep the reactive wavepacket inside 
the box. The forward propagation is terminated when 
the overlap between the trial wavefunction and the prop- 
agated wavefunction is maximized. 

Starting from the trial wavefunction of Fig. ^ the 
pulse shaping scheme was used to obtain the optimal 
wavepacket in the reactant channel leading to desorp- 
tion of the molecule in the fifth vibrational state. Fig. 
shows the accumulated flux of desorbing molecules as 
a function of the iteration number. State-resolving the 
accumulated flux enables to calculate the percentage of 
desorbing molecules in the fourth and fifth vibrational 
states Cf. Fig. |5l A monotonic increase of the yield in 
V — A is observed. 

The optimal reactant wavefunction is shown in Fig. 01 
(left). The wave function of the matter wave is projected 
out for the equilibrium distance between the chemisorbed 
atom and the surface (0.916A). The ID wave packet is 
found to be multi-pulsed reconfirming previous studies, 
that control is obtained by destructive or constructive 
quantum interferences. The momentum representation is 
obtained by performing a Fourier transform of the one- 
dimension matter-wave. The momentum components are 
shown in Fig. 0] Several momentum components can 
be observed which are needed for creating the necessary 
constructive quantum interference in the fifth vibrational 
state and destructive in all the other vibrational states. 

The source of the matter-waves, the BEC condensate, 
has to be positioned very close to a surface. Such a de- 
vice is realized in the so-called atom chips 
or surface micro-traps 16. 17]. In the atom chips for ex- 
ample the BEC has been placed a few hundred microns 
above a metal surfaceTi^. 

The experimental realization of the control scheme de- 
pends on the the ability to shape matter- waves. Four- 
wave mixing I IS'l of matter waves has been demonstrated 
which can lead to different momentum components of 
the matter-waves. A more comprehensive solution could 
follow ideas from optical pulse shaping. A matter-wave 
diffraction grid has been demonstrated^, H^- A setup 



of two diffraction grids and a phase shift device would 
constitute the ultimate matter- wave pulse shaper. Such 
a device would make the suggested control scheme real- 
izable. 
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